Packing and Self-assembly of Truncated Triangular Bipyramids 
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Motivated by breakthroughs in the synthesis of faceted nano- and colloidal particles, as well as 
theoretical and computational studies of their packings, we investigate a family of truncated trian- 
gular bipyramids. We report dense periodic packings with small unit cells that were obtained via 
numerical and analytical optimization. The maximal packing fraction (/)max changes continuously 
with the truncation parameter t. Eight distinct packings are identified based on discontinuities 
in the first and second derivatives of (/)max(t). These packings differ in the number of particles in 
the fundamental domain (unit cell) and the type of contacts between the particles. In particular, 
we report two packings with four particles in the unit cell for which both 0max(t) and (/)max(^) 
are continuous and the discontinuity occurs in the second derivative only. We also perform self- 
assembly simulations for large system sizes. Only one out of eight packings is found to assemble 
in the thermodynamic limit. In addition, the degenerate quasicrystal reported previously for tri- 
angular bipyramids without truncation [Haji-Akbari et a/., Phys. Rev. Lett. 107: 215702 (2011)] 
assembles for truncations as high as 0.45. The self-assembly propensities for the structures formed 
in the thermodynamic limit are explained using the isoperimetric quotient of the particles and the 
coordination number in the disordered fluid and in the assembled structure. 



I. INTRODUCTION 



Hard core interactions are idealized representations of 
excluded volume effects that are prevalent in dense states 
of matter such as liquids and solids, and have been shown 
to approximate short-range repulsions [1]. In atomic and 
molecular systems, such repulsions arise when atoms and 
molecules get close enough for their electronic shells to 
overlap and repel. It is therefore no surprise that hard 
spheres have been successful models of dense liquids and 
solids for decades [2-6]. At larger length scales, hard 
core interactions are used to predict the thermodynamic 
behavior of concentrated colloidal suspensions since col- 
loidal particles can, under certain circumstances, behave 
like hard particles. For example, long-range electrostatic 
repulsions between charged nanoparticles are screened by 
ionic solutions, effectively becoming short-range repul- 
sions [7], which can be approximated with hard core in- 
teractions [1]. The phases expected for such colloidal sys- 
tems can be predicted through theoretical and computa- 
tional studies of the corresponding hard particle system. 
Recent breakthroughs in the synthesis of colloidal parti- 
cles of various shapes [8-20] have spurred an increased 
interest in studies of hard particle systems. 

A remarkable fact about hard particles is their abil- 
ity to spontaneously form ordered phases at sufficiently 
high packing fractions. These phases are stabilized by 
entropy alone and are observed in computational studies 
of various hard particle shapes including spheres [2, 3], 
spherocylinders [21], thin disks [22, 23], ellipsoids [24- 
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28], dumbbells [29], tetrahedra [30, 31], triangular bipyra- 
mids [32], superballs, cubes and octahedra [33, 34], snow- 
man particles [35], squares [36], space-filling polyhe- 
dra [37-39] and many other polyhedra [40] . The prepon- 
derance of optical, electrical, magnetic and mechanical 
properties of ordered structures formed from such parti- 
cles [41, 42] makes knowledge and prediction of their ex- 
pected thermodynamic assemblies of particular current 
interest. 

The feasibility of disorder-order transitions can be un- 
derstood from a packing perspective as all permissible 
configurations of hard particle systems are packings of 
the corresponding shape. More particularly, the struc- 
ture that is thermodynamically stable in the limit of 
infinite pressure is the packing (or, if there are several 
equally dense packings, is among the packings) with the 
largest possible packing fraction. This is seen from the 
Gibbs free energy G(P, T) = PV - ST of the hard par- 
ticle system. As the entropy is bounded from below and 
decreases with increasing pressure, the Gibbs free en- 
ergy is dominated by the PV term as P ^ oo. The 
structure with the smallest volume (largest packing frac- 
tion) will eventually be stable. Solving the packing prob- 
lem is therefore equivalent to identifying the equilibrium 
structure of the corresponding hard particle system in 
the infinite-pressure limit. All known densest packings 
of two- and three-dimensional objects are ordered [43], 
which necessitates a disorder-order transition at some 
finite pressure and density. However, it is also known 
that the ordered structures formed by many hard par- 
ticle systems at intermediate densities can be geomet- 
rically unrelated to their corresponding densest pack- 
ings [31, 32, 40, 44]. 

One hard particle with a rich and unconventional phase 



behavior is the hard tetrahedron. The densest known 
packing of tetrahedra is a double-dimer lattice with four 
particles in its fundamental domain forming two pairs 
or 'dimers' [45-47]. At intermediate densities, however, 
tetrahedra self-assemble into a dodecagonal quasicrys- 
tal [30]. The (3.4.3^.4) approximant of the quasicrystal 
is only slightly less dense than the double-dimer lattice 
and is more stable at intermediate densities [31]. A sim- 
ilar trend is observed for the closely related hard trian- 
gular bipyramid, which is a dimer of tetrahedra. Tri- 
angular bipyramids form a degenerate dodecagonal qua- 
sicrystal, a structure that is identical to the tetrahedron- 
based quasicrystal on the monomer level but random 
in the pairing of constituent tetrahedra into dimers in 
the nearest neighbor network. At intermediate densities, 
both the degenerate quasicrystal and its approximant are 
again more stable than the densest packing of triangular 
bipyramids, which is structurally identical to the double- 
dimer lattice of tetrahedra [32] . 

Damasceno et at. studied the self assembly and densest 
packing behavior of a family of truncated tetrahedra [44] . 
The building blocks were obtained by truncating the ver- 
tices of a regular tetrahedron. The particle geometry 
ranged from a perfect tetrahedron with truncation pa- 
rameter t = to a perfect octahedron with t = 1. It 
was shown that the formation of the quasicrystal for the 
tetrahedron at intermediate packing fractions is robust 
and occurs for truncations up to t < 0.45. Several un- 
expected structures self-assemble at higher truncations, 
including diamond, /3-tin, and a crystal isostructural to 
high-pressure lithium. A previously unreported space- 
filling polyhedron was obtained for t = 1/2. Most no- 
tably, the assembled structures of truncated tetrahedra 
are typically distinct from the densest packings obtained 
in numerical compressions of small unit cells [44] . In ad- 
dition to the potentially complex phase behavior, study- 
ing truncated polyhedra is relevant from a practical per- 
spective. Experimentally synthesized nano- and colloidal 
particles are not perfect, but typically have various types 
of imperfections [48] , one of which can be truncated ver- 
tices [13]. Studying the role of truncation allows us to 
investigate the robustness of the self-assembly process to 
such imperfections. 

Motivated by these considerations, we use techniques 
similar to those used in Refs. [30, 32, 44, 47] to inves- 
tigate the packing and self-assembly of hard truncated 
triangular bipyramids. We study the densest packings 
with small repeating unit cells using Monte Carlo sim- 
ulations and further analytical and numerical optimiza- 
tion. Eight distinct packing types, six with two particles 
in the unit cell and two with four particles in the unit 
cell, are identified for different values of truncation based 
on discontinuities in the first and the second derivatives 
of 0max(^)- The two four-particle packings are observed 
for a narrow range of truncations and are separated by a 
discontinuity in ^max(^)- We also simulate large systems 
of truncated triangular bipyramids and observe that only 
one out of eight packing types forms in simulation. We 



observe that the quasicrystal formation is robust and oc- 
curs for truncations as high as 0.45. Finally, we explain 
the self-assembly potential of different building blocks by 
their isoperimetric quotient and the coordination number 
of their disordered fluid and find that only the structures 
that are locally similar to their disordered fiuid tend to 
form in the self-assembly simulations of the bulk system. 
This paper is organized as follows. Geometrical no- 
tations and definitions are given in Section II A. Tech- 
nical details of Monte Carlo simulations are presented 
in Section II B. Section II C outlines our formulation of 
the constrained packing problem and the approach used 
for its solution. The solutions of the packing problem 
for truncated TBPs are thoroughly discussed in Sec- 
tion III A. Section IIIB discusses the self-assembly simu- 
lations while the role of the local structure and the parti- 
cle shape on self-assembly are discussed in Section IIIC. 
And Section IV is reserved for discussions and concluding 
remarks. 



II. METHODS 

A. Geometric notations and definitions 

Let 7^ be a perfect (non-truncated) triangular bipyra- 
mid (TBP) centered at the origin. It is the convex hull 
of its five vertices 



o = (+2,+2,+2), 
P = (+2,-l,-l), 
q=(-l,+2,-l), 
r=(-l,-l,+2), 
s = (-2,-2,-2). 



(la) 
(lb) 
(Ic) 
(Id) 
(le) 



To truncate the TBP, we replace these vertices by new 
vertices positioned along the edges as specified by a trun- 
cation parameter t. Two new vertices are positioned uni- 
formly and symmetrically on each edge in-between the 
old vertices, t = 0, and the edge mid centers, t = 1. 
The truncated triangular bipyramid (tTBP) Vt is then 
defined as the convex hull of the following 18 vertices 



op = (1- |)o+|p, 
oq= (1- |)o+|q, 
or=(l-|)o+|r. 



po=(l-|)p+|o, (2a) 
qo=(l-|)q+|o, (2b) 
ro = (l-|)r+|o, (2c) 



sp = (1 - 


-§)s+|p, 


ps = (1 - 


-§)p+§s, 


(2d) 


sq= (l- 


-|)s+iq, 


qs = (1 - 


-§)q+|s, 


(2e) 


sr=(l- 


-|)s+|r, 


rs = (1 - 


-|)r+|s, 


(2f) 


qr = (1 - 


-|)q+|r, 


rq = (1 - 


-|)r+|q, 


(2g) 


rp = (1 - 


-|)r+|p, 


pr = (1 - 


-|)p+|r, 


(2h) 


pq= (l- 


■|)p + |q> 


qp = (1 - 


-|)q+|p. 


(2i) 



We denote a vertex of Vt by V^[x], where x is any of 
the points given in Eq. (2). An edge is denoted either 




^ = 2/3 



t = l 



FIG. 1: (a) A truncated triangular bipyramid consists of two 
regular tetrahedr a joined together at one face, truncated sym- 
metrically at the five vertices. Polar triangles, equatorial tri- 
angles, and peripheral hexagons are denoted by A, B, and C, 
respectively, (b) Truncated triangular bipyramids for differ- 
ent values of the truncation parameter t. 



by its endpoints or by any two points that lie along the 
line segment connecting its endpoints. Similarly, a face 
is denoted by any three points that are coplanar with the 
face. For instance, an edge with the endpoints op and po 
is denoted by E'[o,p] and the hexagonal face containing 
the points op, po, pq, qp, qo, and oq by F[o, p, q]. 

Every tTBP has three types of faces that are distin- 
guished in Fig. la. The two triangular faces perpendic- 
ular to the three-fold axis of the particle are referred to 
as polar triangles and are labeled A in the figure. The 
six triangles with an edge parallel to the three-fold axis 
are called equatorial triangles and are labeled B. The re- 
maining six hexagonal faces are referred to as peripheral 
hexagons and are labelled C. 

The volume and the surface area of Vt are given by 



h') 



V{t) = V, (1 

Sit) = Sp \l-U6-V6)i' 



(3) 

(4) 



where T^ = 18 and Sp = 27>/3 are the volume and the 
surface area of 'P, respectively. Fig. lb depicts tTBPs for 
a few special values of t. 



B. Monte Carlo simulations 

We perform Monte Carlo (MC) simulations of tTBPs 
with truncations ranging from zero to unity. The edge 
length of a perfect non-truncated TBP is chosen as 



the length scale and pressure is measured as P* = 
PV{t)/kBT. Overlaps between tTBPs are detected us- 
ing the Gilbert- Johnson-Keerthi (GJK) algorithm [49]. 
Each MC cycle is comprised of one trial move per parti- 
cle, which is either a trial translation or a trial rotation 
with equal probabilities. All simulations are performed 
in boxes with periodic boundary conditions. 

Two types of MC simulations are carried out in this 
study. In order to identify dense small-cell packings of 
tTBPs, isotension MC simulations are performed in small 
unit cells, i.e. in boxes containing up to eight particles. 
Such simulations include an additional trial move per cy- 
cle wherein both the size and the shape of the simulation 
box are varied. The box move is accepted or rejected 
according to the conventional Metropolis criterion in the 
NPT ensemble. All simulations with small boxes start 
at a low pressure of P* = 1. The pressure is then expo- 
nentially increased every 50 000 MC cycles in a total of 
2 X 10^ cycles to a final value of P* = 10^. For each trun- 
cation and system size, 100 independent simulations are 
performed and the densest packing within that 100-run 
sample is recorded. The maximum packing fractions ob- 
tained from these 100 simulations are within less than 0.5 
per cent of the Umax's obtained from analytical and/or 
numerical solution of the packing problem. 

In order to assess the self-assembly potential of differ- 
ent tTBPs, MC simulations are performed in the canon- 
ical ensemble for packing fractions ranging from 55% to 
69%. We start each simulation from a configuration of 
N = 2 048 tTBPS rapidly compressed to a packing frac- 
tion of ~ 70% and rescaled to the appropriate packing 
fractions. These simulations are performed for a mini- 
mum of 5 X 10^ MC cycles. 



C. The packing problem 



A collection of shapes {7^^}^^ C M'^ is called a pack- 
ing if their interiors are disjoint. For every packing, the 
packing fraction (or packing density) is defined as 



lim sup 



1 



Vol[P(r)] 



^Vol[7^,nP(r)], (5) 



with B{r) = {a; G M'' : ||a;|| < r}. For a periodic packing 
of congruent objects, Eq. (5) takes the form 



NVn 
Vb ' 



(6) 



where N is the number of particles in the unit cell. Vjz 
and Vb are the volumes of the individual particle and the 
unit cell, respectively. 

The packing problem searches for the densest possi- 
ble packing and can be stated as follows: What is the 
arrangement of the shapes {7li}i that maximizes (/) as 
defined in Eq. (5)? In its most general formulation, with- 
out a confining box, the packing problem is an optimiza- 
tion problem with an infinite number of variables, and 



thus difficult to solve. However, since packings in spa- 
tially restricted systems depend only on a finite number 
of parameters, the local maxima of Eq. (5) can often be 
determined numerically or (in a few cases) even analyti- 
cally. For instance, the densest lattice packings of three- 
dimensional compact convex shapes have been obtained 
using the classical method of Minkowski [50] for regular 
tetrahedra [51] and regular octahedra [50]. Betke and 
Henk developed an efficient computer algorithm for the 
determination of dense lattice packings of 3-polytopes, 
and applied it to Platonic and Archimedean solids [52]. 
Local maxima have been obtained for periodic non-lattice 
packings with a few particles in the unit cell for several 
three-dimensional objects [43, 45, 47, 53, 54]. 

In this work, we use the following procedure proposed 
by Chen et al. [47] . Given a dense packing obtained from 
MC simulation, we identify all contacts between the par- 
ticles in the unit cell. These contacts are then expressed 
as intersection equations, which act as constraints for the 
suboptimal optimization problem: 



minimize V{qi,q2,--- .Qn) 
subject to 4((7i,^2,--- ,Qn) 



0, k = 1,2, 



, m. 



(7) 
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FIG. 2: Packing fraction (/)max(t) of the densest packings 
and self-assembled arrangements of tTBPs. The upper sym- 
bols are the maximum packing densities obtained from MC 
simulations, and the curves are calculated from the analyti- 
cal expressions given in Appendix A or numerical solutions 
outlined in Appendix B. In addition, the lower symbols cor- 
respond to the lowest densities at which each phase forms in 
self-assembly simulations. 



where V{qi^q2r - ■ ^Qn) ^^ the volume of the periodic box. 
The Qi correspond to the n degrees of freedom of the 
packing (examples are lattice vectors, positions and ori- 
entations of particles in the unit cell) and Ik = are 
the m intersection constraints that need to be satisfied 
by the packing. Eq. (7) is solved analytically or numer- 
ically to obtain the densest packing consistent with the 
constraints. 



III. RESULTS 



TABLE I: The eight types of densest tTBP packings. We list 
the numbers of independent intersection constraints for each 
packing. 



Type 


P2I 


P22 


P23 


P4I 


P42 


P24 


P25 


P26 


Vertex-face 


1 








3 


5 











Edge-edge 


2 


2 


2 


2 


2 


4 


2 


3 


Edge- face 











2 


2 











Face- face 


8 


8 


8 


3 


3 


5 


8 


6 


Total 


11 


10 


10 


10 


12 


9 


10 


9 



A. Dense packings 

The packing fractions determined from simulations 
with small boxes are reported in Fig. 2. As with any con- 
tinuous deformation, 0max(O changes continuously with 
the truncation t. We identify eight distinct packing types 
by noting six discontinuities in the first derivative ^max(^) 
and one discontinuity in the second derivative ^max(^)- 
All packings have densities that exceed 84%, which shows 
that tTBPs are generally efficient packers. The tTBP 
with truncation t = | is a space-filling polyhedron pre- 
viously indexed as 14-III by Goldberg [55]. Eq. (7) was 
formulated for all the packings. Details can be found in 
the Appendices A and B. We label packings as P^^, where 
n specifies the number of particles in the unit cell and i is 
a running index characteristic of the packing type. The 
packings and their interaction constraints are listed in 
Table I. 



1. Two-particle packings 

Six of the eight densest packings have two particles in 
the minimal unit cell. Eq. (7) was solved analytically 
for all six packings and the full mathematical description 
of each packing is given in Appendix A. In each packing, 
the orientations of the two tTBPs are related by inversion 
symmetry, which means they are Kuperberg pairs [56]. 
The layers formed by each particle will be denoted by 
+ and - for the original orientation and the inversion, 
respectively. 

For small truncations, we find a packing that is similar 
to the TBP crystal proposed in [47] except for a gradual 
shear to optimize the packing of slightly truncated TBPs. 
This packing type, which we denote by P2I (Fig. 3a), 
has eleven independent intersection constraints that are 
given in Eq. (A2), two more than the original TBP crys- 
tal [47]. The face-to- face contacts in the P2I packing are 
all among peripheral hexagons and neither the polar nor 




FIG. 3: The packing P2I for the truncation t = |. (a) The 
polar triangle F[op,oq, or] of particle B touches the pe- 
ripheral hexagon F[o,p,q] of particle A through the edge 
E[oq_,or]. (b) The central blue particle has face-to-face con- 
tacts with eight of its neighbors, all of which are in an opposite 
layer, (c) Particles with identical orientations touch through 
edge-to-edge contacts only. 




FIG. 4: The packing P22 for the truncation t = ^ (a) and 
the packing P23 for the truncation t = | (b). (a) In P22, the 
faces F[op,oq, or] of the particles A and B touch and form 
a star-like pattern. The axial edge £^[ro, rs] of particle A 
and the peripheral edge ^[pq, qp] of particle C touch, (b) In 
P23, polar faces of particles A and B touch forming a star-like 
pattern. No axial edge of P23 touches an edge of a neighboring 
particle. 



the equatorial triangles touch any other faces. All face- 
to-face contacts are between particles in opposite layers 
(Fig. 3b) and particles with identical orientations touch 
through edge-to-edge contacts only (Fig. 3c). The list 
of face-to- face contacts in P2 1 is identical to that in the 
TBP crystal. The maximum packing fraction for P2I is 
given by 



lit) = 



500(8 - 3^3) 
4671 - 30t - 25t2 ■ 



(8) 



The next two packing types, P22 (Fig. 4a) and P23 
(Fig. 4b), are characterized by star-like patterns formed 
by the face-to-face touching of polar triangles of parti- 
cles in opposite layers (Fig. 4). This is in contrast to 
P2I where the polar triangles only touch the peripheral 
hexagons through one of their edges (Fig. 3a). Both 
P22 and P23 have a total of ten independent intersec- 
tion constraints given in Eq. (A4) and Eq. (A6), respec- 
tively. The major difference between these two otherwise 
similar packings is the existence of an edge-to-edge con- 
tact in P22, where an axial edge (an edge parallel to the 
three-fold axis of the particle) touches a peripheral edge 
(Fig. 4a). The maximum packing fractions of P22 and 
P23 are given by 



2,2 (t) = 


283 + 131t-319t2+64t3' 
27(8 - 3^3) 


(9) 
(10) 


<p2Ait) = 

<i>2At) = 


2(8 - 3^3) 


2,3(t) = 


3(2-i)(4 + 4t-7t2)' 


2(l+t)(119-80t + 8t2)- 


28(8 - 3^3) 



The next two types of two-particle packings, P24 
(Fig. 5a-b) and P25 (Fig. 5c-d), are characterized by the 



contact of polar triangles of identically oriented particles 
(Fig. 5a,c). This is unlike P22 and P23 where star-like 
patterns form between polar triangle of particles in al- 
ternate layers (Fig. 4). In both P24 and P25, particles 
arrange into a hexagonal lattice in each layer with their 
equatorial triangles touching. Face-to-face contacts be- 
tween particles in opposite layers occur through periph- 
eral hexagons. Therefore, the + and - layers are shifted 
with respect to one another and the entire structure 
is isostructural to a stretched hexagonally close-packed 
(hep) lattice (Fig. 5b, d). 

There are two major structural differences between 
these two packings. The face-to- face contacts between 
polar triangles are perfect in P25 (Fig. 5c), while the 
touching polar triangles are shifted in P24 (Fig. 5a). The 
contacts between equatorial triangles are also different. 
In P24, equatorial triangles of triplets of particles come 
together, like for example B, C and D in Fig. 5b. Yet, 
in P25 one of each of the three particles in those triplets 
start distancing from the other two as the truncation pa- 
rameter increases. This leads to an elongation of the 
intra-layer hexagonal lattice along one of its lattice vec- 
tors (Fig. 5d). Overall, P24 and P25 have nine and ten 
independent intersection constraints that are given by 
Eq. (A8) and Eq. (AlO) respectively. The maximum 
packing fractions are given by 

(11) 

(12) 

The boundary between the P24 region and the P25 re- 




FIG. 5: The packing P24 for the truncation t = | (a,b) and 
the packing P25 for the truncation t = | (c,d). (a,c) The side 
view reveals that the face F[op, oq, or] of particle A touches 
the face i^[sp, sq, sr] of particle B. The contact is shifted in (a) 
but not in (c). Equatorial triangles of particles B, C and D 
come together in (b,d). It is visible in the top view (d) that 
particle D shifts away from the (B,C,D) triplet as t becomes 
larger. 



gion at t = I corresponds to a tTBP that tiles Euclidean 
space (Fig. 10a). This space-filling polyhedron is known 
and was indexed as 14-III by Goldberg in 1979 [55]. We 
rediscover this polyhedron in our study and observe that 
the corresponding tiling forms in self-assembly simula- 
tions of the same building block (see below in Fig. lOb-c). 

The last type of two-particle packings, P26, is charac- 
terized by the absence of contacts between two consecu- 
tive -\- and - layers. Like P22 and P23, polar triangles 
of particles in opposite layers intersect, but the contact 
is less perfect and the particles tend to have a relatively 
large lateral shift (Fig. 6b). This leads to a staircase 
arrangement for each layer. Like P24, the equatorial 
triangles of particles in identical layers are face-to-face 
(Fig. 6a). The peripheral hexagons also touch the cor- 
responding hexagons from opposite layers. This packing 
has a total of nine independent intersection constraints 




FIG. 6: (a) Top and (b) side view of the packing P26 for the 
truncation t = || . 



given by Eq. (A12). Its packing fraction is given by 

20(8 - 3t^) 



/^2,6 



it) = 



3(2-i)(52 + 20t-35t2)' 



(13) 



Further information about these packings, including 
analytical expressions for lattice vectors and intersection 
constraints, are given in Appendix A. 



2. Four-particle packings 

The two packing types observed for the intermediate 
truncations ts ^ 0.5010 < t < t^ ^ 0.5776 have four 
particles in the fundamental domain. The basic building 
block of these packings is a (non-convex) dimer of tTBPs 
with partially touching peripheral hexagons (Fig. 7a). 
Each dimer is a Kuperberg pair as the orientations of 
its constituent tTBPs are related by inversion. The fun- 
damental domain contains two such dimers that are ro- 
tated with respect to one another at an angle of about 
60 degrees. All face-to- face contacts are between parti- 
cles in Kuperberg pairs with identical orientations: Po- 
lar triangles of particles with identical orientations touch 
(Fig. 7c), while four out of the six peripheral hexagons of 
each particle touches the peripheral hexagons of neigh- 
bors with opposite orientations, i.e. the neighbors with 
orientations related to the central particle by inversion 
(Fig. 7d). In contrast, the contacts between particles in 
separate dimer types are either edge-to-edge or vertex-to- 
face (Fig. 7b). There are a total of ten and twelve inde- 
pendent intersection constraints for P4I and P42, respec- 
tively, alongside eight additional symmetry constraints 
given by Eq. (B2). Due to the non-linear nature of the 
intersection constraints, no analytical solution can be ob- 
tained for Eq. (7). Numerical solutions were however, 
obtained and the results are depicted in Figs. 2 and 8. 

These two packings are similar to the SM2 packing 
of hard ellipsoids [26, 27], which has a monoclinic unit 
cell with two ellipsoids forming an angle between them. 
The SM2 phase is denser than the stretched face-centered 
cubic (fee) phase, which is obtained from an affine trans- 
formation of the sphere fee packing [57]. For sufficiently 
large aspect ratios, the SM2 packing is always more sta- 
ble than the stretched fee phase [28]. 
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FIG. 7: Dense packings of tTBPs with four particles in 
the unit cell for the truncation t = ^. (a) The dimer is a 
Kuperberg pair and the basic building blocks of the packings 
P4I and P42. The red and blue tTBPs belong to the same 
dimer, while the green and yellow particles belong to the other 
dimer rotated by an angle of about 60° . The crystal structure 
is depicted from different perspectives in (b-d). 



A remarkable observation about the two four-particle 
packings is that all intersection constraints satisfied by 
P4I are also satisfied by P42. In addition, P42 satisfies 
two additional constraints that correspond to two inter- 
dimer vertex-to-face contacts. This makes the set of all 
feasible packings for P42 a subset of all feasible packings 
for P4I, and leads to a unique transition of the packing 
behavior in the truncation space. While both ^max(^) 
and 0max(^) change continuously from P4I to P42, a dis- 
continuity occurs only in the second derivative (Fig. 8). 
We therefore characterize this transition as a 'second- 
order' transition in the truncation space. Analogously, 
the usual behavior of a discontinuity in ^max(^) would be 
called a 'first-order' transition. 



B. Self-assembly simulations 



FIG. 8: The packing fraction (/)max(t) and its derivative 
0max(^) are continuous at the transition from P4I to P42. 
Magenta squares are the densest packings obtained from 
small-box MC simulations. The blue and the red curves are 
the packing fractions calculated from numerical solution of 
Eq. (7). These curves are depicted in solid lines wherever 
they correspond to the densest packing and dashed other- 
wise. As can be seen in the inset, there is a discontinuity in 
the slope of (/)max(^) at t4 ~ 0.5321, which corresponds to a 
second-order transition in the truncation space. 



TABLE II: Formation of the degenerate quasicrystal of 
tTBPs. (j)min corresponds to the lowest packing fraction at 
which the quasicrystal assembles from the disordered fluid. 
0max is to the maximum packing fraction obtained from a 
numerical compression of the assembled bulk quasicrystal. 



Truncation 


(/>min 


MC Cycles 


0niax 


0.1 


57% 


4x 10^ 


81.21% 


0.2 


57% 


6 X 10^ 


80.93% 


0.3 


58% 


9 X 10^ 


79.74% 


0.4 


60% 


15 X 10^ 


81.31% 


0.45 


61% 


28 X 10^ 


80.06% 



The tTBP family assembles far fewer ordered struc- 
tures than the truncated tetrahedron family studied in 
Ref. [44] . We observe the assembly of only two distinct 
crystal structures. A dodecagonal quasicrystal forms for 
truncations t < 0.45. We already reported the forma- 
tion of this quasicrystal in a system of non-truncated 
TBPs [32]. At higher truncation, P25 forms in the 
range 0.66 •• • < t < 0.80. For the remaining trun- 
cations, no spontaneous ordering was observed on the 
(long) timescale of our simulations. 

We analyze the structures obtained from our simula- 
tions using the tools outlined in Ref. [32]. Table. II lists 
the minimum packing fraction at which the degenerate 
quasicrystal forms for each truncation as well as the num- 
ber of MC cycles it took for it to form in our simulations. 
As expected and as also observed with imperfect tetra- 
hedra [48], it becomes increasingly difficult for tTBPs 
to form the quasicrystal as the truncation increases and 



higher packing fractions and longer simulation times are 
necessary for crystallization. 

To compare the structural quality of the quasicrystals 
to the degenerate quasicrystal obtained for zero trunca- 
tion, we calculate the diffraction patterns for the assem- 
bled structures compressed to their maximum packing 
fractions (so-called "inherent structures" [58]). Due to 
the degenerate nature of the quasicrystal, the scatterers 
are placed not at the centroids of TBPs, but instead at 
the centroids of each of the two tetrahedra that comprise 
each non-truncated TBP. The resulting diffraction pat- 
terns resemble the ones obtained for a system of hard 
tetrahedra [30]. The twelvefold symmetry in the diffrac- 
tion pattern, which is the signature property of a do- 
decagonal quasicrystal, is readily visible in all diffraction 
patterns (Fig. 9). It is notable that neither the qual- 
ity of the diffraction patterns nor the maximum (post- 
compression) packing fractions of the quasicrystal(s) are 



significantly different from the degenerate quasicrystals 
formed by untruncated TBPs. This suggests that the 
structural quality of the quasicrystal is only slightly com- 
promised as a result of truncating the building blocks, 
even though it gets gradually harder to assemble. 

For truncations 0.66 •• • < t < 0.80, the system as- 
sembles into a simple crystalline structure identical to 
the P25 packing (described in section IIIAl). Fig. 10b 
shows the final snapshot of an isochoric simulation of 
2 048 tTBPs with t = § at ^ = 62%. In Fig. 10a, we 
show the analytically constructed P25 packing for the 
same truncation, which is one of Goldberg's space-filling 
tilings [55]. In both panels, the location of the parti- 
cles in the + and - layers are highlighted with yellow 
and cyan lines, respectively. Fig. 10c shows the diffrac- 
tion pattern calculated when viewed along the vector per- 
pendicular to the layers; the scatterers are positioned at 
the cetroids of tTBPs and not the constituent truncated 
tetrahedra. The observed sixfold symmetry is consistent 
with the elongated hep structure of the P25 packing. For 
comparison, the truncated tetrahedron system at t = | 
assembles into a diamond crystal [44]. The packing ob- 
served in this study is not degenerate to the diamond 
crystal observed for truncated tetrahedra. 



C. The role of particle shape and local order 

Except for P25, none of the dense packings found in 
this study form in self-assembly simulations. As men- 
tioned in Section I, the densest packing of every hard par- 
ticle will be thermodynamically stable eventually for suf- 
ficiently large pressures. As previously shown for tetra- 
hedra [31], truncated tetrahedra [44], triangular bipyra- 
mids [32], and several other shapes [40], hard particles 
can assemble at intermediate densities into structures 
that are very different from their densest packings. While 
it has been observed that the assembled packings are 
typically similar to the dense disordered fluid on a lo- 
cal level [40] , no such relation is known between the fluid 
or the assembled packings and the densest packings. 

There are several ways to quantify the similarity of two 
packings on a local level. The simplest approach is to use 
the coordination number (CN) by integrating the radial 
distribution function ^(r). 



CN 



47rpnr'^g{r)dr. 



(14) 



with a cut-off tq at the minimum between the first and 
second peak. Here pn is the average number density. It 
has been demonstrated recently by Damasceno et al. [40] 
that there is a nearly perfect correlation between the 
CN of the dense disordered fluid and the CN of the self- 
assembled structure. We calculate CN for the dense dis- 
ordered fluid, the densest packing, and the self- assembled 
structures at ^ = 60%. The results are presented in Ta- 
ble III. All assembled structures have CNs that are close 
to that of the disordered fluid at the same density. It is 



noteworthy that the proximity of CNs of the fluid and 
a particular packing is not a sufficient condition for self- 
assembly to occur since CN is only a spatially averaged 
measure of local order and does not take into account 
the rotational anisotropy of the nearest-neighbor shell. 
Moreover, CN can be problematic for certain elongated 
shapes, such as ellipsoids since it can underestimate the 
local coordination shell. However when there are stark 
differences in the CNs of the fluid and the densest pack- 
ing, for example at t > 0.85, the corresponding densest 
packing never forms in our self-assembly simulations. 

Further, we observe an abrupt change in the CN of 
the disordered fluid at around t = 0.6. This change can 
be understood by directly inspecting the behavior of the 
radial distribution function of the disordered fluid. As 
can be seen in Fig. 11, the first peak of g{r) grows and 
broadens for truncations from t = 0.66 to t = 0.80, which 
increases CN for these truncations. Near t = 0.85 a small 
second peak appears at r ^ 0.8, which reduces CN by our 
definition and may explain why P25 fails to form for trun- 
cations beyond 0.8, despite being the densest packing for 
truncations as high as 0.95. A similar peak-splitting oc- 
curs in P25 near t = 0.90 (Fig. 12) which leads to an 
abrupt decrease in CNpacking from 5.1 to 2.7. This is in 
line with our earlier observation that layers are stretched 
along one lattice direction as t increases. The shifting 
of particles in the triplets of tTBPs sharing peripheral 
hexagons (explained in Section IIIAl) makes the struc- 
ture less uniform, which might also be a contributing 
factor in its kinetic inaccessibility. 

Damasceno et al. [40] also find a correlation be- 
tween the isoperimetric quotient (IQ) of a crystal- forming 
building block and the type of crystal that it assembles. 
Crystal formers with intermediate IQs tend to form non- 
Bravais lattices while those with large IQs are more likely 
to form face-centered and body-centered cubic lattices. 
This is consistent with our findings in Table HI as the 
particles forming the quasicrystal all have mostly inter- 
mediate isoperimetric quotients and the particles forming 
the P25 packing have IQs compatible with simple crystal 
formation. 



IV. DISCUSSION AND CONCLUSIONS 

We found that the formation of the degenerate do- 
decagonal quasicrystal from truncated triangular bipyra- 
mids is robust and can be observed for truncations as 
high as 0.45. This is in line with earlier studies of the 
truncated tetrahedron system, and might be a conse- 
quence of the close local structural similarity of the qua- 
sicrystal and the disordered fluid. The quasicrystalline 
phase has been shown to be more favorable thermody- 
namically than the densest known packings for tetrahe- 
dra [31] and triangular bipyramids [32] at intermediate 
densities. Whether the same is true for tTBPs or trun- 
cated tetrahedra, is beyond the scope of this work and 
can be the subject of future studies. 
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FIG. 9: Diffraction images of the degenerate dodecagonal quasicrystals of tTBPs for the truncations (a) t = 0.1, (b) t 
(c) t = 0.3, (d) t = 0.4, and (e) t = 0.45. The sharpness and intensity of the peaks decreases slowly with truncation. 



0.2, 



TABLE III: Coordination number (CN) for the dense disordered fluid, the self- assembled structure, and the densest packing at 
different values of truncation. The isoperimetric quotient (IQ) of a particle p is defined as IQ = SGnVp /Sp. It is a measure of 
the particle sphericity [40]. 



Truncation 


IQ 


Densest Packing 


Self-assembled Structure 


CNfl^id 


CNself-assembled 


CJN packing 


0.10 


0.3623 


P2I 


Quasicrystal 


1.9 


1.9 


1.7 


0.20 


0.3735 


P2I 


Quasicrystal 


1.9 


1.9 


1.5 


0.30 


0.3913 


P22 


Quasicrystal 


2.0 


1.9 


4.0 


0.40 


0.4150 


P23 


Quasicrystal 


2.0 


1.9 


2.0 


0.45 


0.4290 


P23 


Quasicrystal 


2.0 


1.9 


1.9 


0.50 


0.4444 


P23 


- 


2.0 


- 


1.9 


0.55 


0.4611 


P42 


- 


2.0 


- 


1.6 


0.60 


0.4791 


P24 


- 


2.2 


- 


3.9 


0.67 


0.5061 


P25 


P25 


4.6 


3.8 


3.8 


0.70 


0.5182 


P25 


P25 


4.5 


3.9 


3.9 


0.80 


0.5600 


P25 


P25 


4.0 


4.2 


4.2 


0.85 


0.5808 


P25 


- 


1.8 


- 


5.1 


0.90 


0.6004 


P25 


- 


1.7 


- 


2.7 


0.95 


0.6176 


P25 


- 


1.5 


- 


2.5 


1.00 


0.6304 


P26 


- 


1.4 


- 


2.6 



The question of the kinetics of quasicrystal formation 
in these systems is undoubtedly interesting and open for 
future exploration. Our studies suggest that the self- 
assembly process is robust and is not affected signifi- 
cantly by small symmetric truncations of the vertices of 
the building blocks. This robustness is especially promis- 
ing from an experimental viewpoint, because it is diffi- 
cult to produce geometrically perfect shapes on nanome- 
ter to micron scales. Our finding hence increases the 
prospects of observing a colloidal quasicrystal for imper- 
fectly shaped triangular bipyramids. 

In all earlier explorations of densest packings obtained 
for a continuously transforming shape [33, 44, 54, 59, 60], 
packing types adjacent in the density plot are conjoined 
via discontinuities in the first derivative of ^max with 
respect to the variable used for the continuous trans- 
formation of that shape. All two-particle packings in 
this study are indeed characterized by discontinuities in 
</^max(0- I^ contrast, the transition from P4I to P42 is 
unique in the sense that a discontinuity occurs in ^^axl^)' 
but not in (/>max(^)- With breakthroughs in the synthe- 



sis of faceted nanoparticles that has spurred interest in 
studying packings for continuously transforming build- 
ing blocks [33, 61], it is likely that more second- and 
higher order transformations will be discovered for other 
families of continuously transforming particle shapes. It 
is noteworthy that these transformations- whether they 
are first-order or second-order- can be experimentally re- 
alized in systems of reconfigurable particles that change 
geometry as a result of external stimuli such as change 
in ionic strength, pH, temperature, or pressure [62-65]. 

Except for the degenerate quasicrystal, none of the 
ordered phases formed in self-assembly simulations of 
tTBPs is degenerate to a structure formed by the trun- 
cated tetrahedron system. This is not surprising since a 
tTBP is not the union of two truncated tetrahedra with 
the same truncation, and the bulkiness around its equato- 
rial triangles makes it geometrically distinct. The added 
bulkiness is relatively small for small truncations, which 
explains why the degenerate quasicrystal still forms in the 
tTBP system. For higher truncations, however, a tTBP 
is so different from a corresponding pair of truncated 
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FIG. 10: (a) Geometrically constructed space- filling packing 
of the tTBP system for t = |. (b) The same structure is 
formed in a simulation of 2 048 tTBPs in an isochoric sim- 
ulation ai (f) — 62%. Colors are chosen randomly, (c) The 
diffraction pattern of (b) calculated from placing scatterers 
at the centers of tTBPs. Little peak broadening and diffuse 
scattering is present, which suggests the assembly is essen- 
tially perfect. 



tetrahedra that none of the structures formed by trun- 
cated tetrahedra is feasible. Truncated TBPs also self- 
assemble into fewer ordered structures than truncated 
tetrahedra, a fact that may be ascribed to the higher 
geometrical anisotropy of tTBPs. 

We explain the observation that only one of the eight 
packing types self- assembles from the disordered fluid 
by using the predictive framework proposed by Dama- 
sceno et al. [40]. Structural differences in the first co- 
ordination shell between the tTBP packings and their 
disordered fluids makes most of these packings less likely 
to assemble from the fluid than a structure with similar 
local packing to that of the fluid, since the latter transfor- 
mation will involve lower energy barriers. A similar phe- 
nomenon is observed for ellipsoids. Oblate and prolate el- 
lipsoids with intermediate aspect ratios have asphericities 
similar to those of the truncated tTBPs considered in this 
study, and- like our packings- are notoriously hard to as- 
semble into their stable SM2 phase from the fluid [26, 27]. 
Further work is needed to better understand the role of 
shape sphericity on self-assembly propensity. 
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FIG. 11: Radial distribution functions of disordered fluids of 
tTBPs with various truncations at packing fraction = 60%. 
The curves are shifted vertically for clarity. 
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Appendix A: Packings of truncated triangular 
bipyramids Avith tAvo particles per unit cell 
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FIG. 12: Radial distribution functions calculated from MC 
simulations of a system of 4 x 4 x 4 unit cells of the P25 packing 
at = 60%. 



In the Appendix, we present analytical expressions for 
the intersection constraints, lattice vectors, and the pack- 
ing fraction of densest packings. In this first section each 
unit cell contains two tTBPs that are related by inver- 
sion. We denote those particles with X and —X, respec- 
tively. The lattice vectors are denoted by a, b and c. 
The offset vector connecting the centroids of X and —X 
is denoted by d. As a result, the centroids of X are given 
by UaS. + n^b + UcC with n^ + ^c + nc = mod 2, and the 
centroids of —X are given by n^a + n^b + ricC + d with 
ria -\- Tie -\- Tie = I mod 2. The volume of the unit cell is 
given by 



V 



det[ b + c, c + a, a + b]. 



(Al) 
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The packing parameters, as denoted by g^ in Section II C, 
are therefore the lattice vectors a, b and c and the off- 
set d. We formulate the intersection constraints for each 
packing and provide lattice and offset vectors for the so- 
lution of Eq. (7) consistent with these constraints. The 
maximum packing fractions are given in the main text 
and thus are not repeated here. The notation used in the 
present work is similar to the notation used in Ref. [47]. 



1. The packing P2I 

There are a total of eleven intersection equations: 

E[o, q] n (^[s, r] + a + b) / 0, (A2a) 

E[o, r] n (^[s, p] + b + c) / 0, (A2b) 

y [rs] n (F[o, p, q] + c - a) / 0, (A2c) 

F[o, p, q] n (d - F[o, p, q] + a) / 0, (A2d) 

F[o, q, r] n (d - F[o, q, r] + b) / 0, (A2e) 

F[o, r, p] n (d - F[o, r, p] + c) / 0, (A2f) 

F[s, q, r] n (d - F[s, q, r] - a) / 0, (A2g) 

F[s, r, p] n (d - F[s, r, p] - b) / 0, (A2h) 

F[s, p, q] n (d - F[s, p, q] - c) / 0, (A2i) 

F[o, r, p] n (d - F[o, r, p] + a + b + c) / 0, (A2j) 

F[s, p, q] n (d - F[s, p, q] - a - b - c) / 0. (A2k) 

The lattice vectors and the offset are given by 



a = 3^(870 + lot, 321 - 25t, -21 + 45t), (A3a) 

b = 3^(-102 - lot, 831 + 25t, 405 - 45t), (A3b) 

c = 3^(282 + 150t, -249 + 145t, 741 + 35t), (A3c) 

d = 3^(38 + lot, 5 + 35t, -25 - 15t). (A3d) 



The maximum packing fraction is given by Eq. (8). This 
packing is valid for < t < ti ~ 0.2544. ti is a root of the 
cubic equation 32000tf - 158825t? + 66310ti - 7117 = 0. 



2. The packing P22 

There are a total of ten intersection equations: 



£;[o,r]n(£;[s,p] + b + c)/0 

E[yo, rs] n {E[p, q] + c - a) / 

F[o, p, q] n (d - F[o, p, q] + a) / 

F[o, q, r] n (d - F[o, q, r] + b) / 

F[o, r, p] n (d - F[o, r, p] + c) / 

F[s, q, r] n (d - F[s, q, r] - a) / 

F[s, r, p] n (d - F[s, r, p] - b) / 

F[s,p,q]n(d-F[s,p,q]-c)/0 

F[op, oq, or] n (d - F[op, oq, or] + a + b + c) / 

F[sp, sq, sr] n (d — i^[sp, sq, sr] — a — b — c) /: 

The lattice vectors and the offset are given by 



36 V 



19t,67-92t,l-23t), 



(A4a) 
(A4b) 
(A4c) 
(A4d) 
(A4e) 
(A4f) 
(A4g) 
(A4h) 
(A4i) 
(A4j) 



(A5a) 



b = ^(16 - lOlt, 103 - 32t, 37 + t), (A5b) 

c = ^(40 + 19t, -17 + 16t, 97 - 23t), (A5c) 

d = ^(-10 + 35t, 5 - lot, 5 - 25t). (A5d) 

The maximum packing fraction is given by Eq. (9). This 
packing is valid for 0.2544 ^ ti < t < t2 = ^. 

3. The packing P23 

There are a total of ten intersection equations: 



£;[o,p]n(£;[s,q]+c + a) /( 

^[o,q]n(£;[s,r]+a + b)/( 

F[o, p, q] n (d - F[o, p, q] + a) / ( 

F[o, q, r] n (d - F[o, q, r] + b) / ( 

F[o, r, p] n (d - F[o, r, p] + c) / ( 

F[s, q, r] n (d - F[s, q, r] - a) / ( 

F[s,r,p]n(d-F[s,r,p]-b)/( 

F[s,p,q]n(d-F[s,p,q]-c)/( 

F[op, oq, or] n (d — F[op, oq, or] + a + b + c) / ( 

F[sp, sq, sr] n (d — i^[sp, sq, sr] — a — b — c) /: ( 

The lattice vectors and the offset are given by 



a 
b 
c 
d 



^(97-14t,40-8t,-17-2t), 
^(l-38t,88 + 4t,67-62t), 
^(37-2t,16-68t,103-26t), 

lot). 



36 V 



36 



- 2t, 16 - 68t, 103 
10t,-10 + 20t,5- 



(A6a) 
(A6b) 
(A6c) 
(A6d) 
(A6e) 
(A6f) 
(A6g) 
(A6h) 
(A6i) 
(A6j) 



(A7a) 
(A7b) 
(A7c) 
(A7d) 



The maximum packing fraction is given by Eq. (10). This 
packing is only valid for | = t2 < t < ta ~ 0.5010. 
There is no analytical formula for ta because there is no 
analytical formula for ^4,1 (t). 



4. The packing P24 



There are a total of nine intersection equations: 



^[op,oq] n (£;[sr, sp] 
£;[qo, qr] n (£^[ps, pr] ■ 
^[ro,rp] n (£;[qs,qp] 
£;[po,pq] n (£;[rs,rq] ■ 



- c + a) / 0, (A8a) 

- b - c) / 0, (A8b) 

- c - a) / 0, (A8c) 

- a - b) / 0, (A8d) 
F[o, p, q] n (d - F[o, p, q] + a) / 0, (A8e) 

F[o, q, r] n (d - F[o, q, r] + b) / 0, (A8f) 

F[o, r, p] n (d - F[o, r, p] + c) / 0, (A8g) 

F[s, p, q] n (d - F[s, p, q] - c) / 0, (A8h) 

F[s, q, r] n (d - F[s, q, r] - a + b - c) / 0. (A8i) 

The lattice vector and the offset are given by 

a = 1(12 - lot, 24 - 22t, 6 - 7t), (A9a) 

b = |(-4t, 24 - 22t, 18 - 13t), (A9b) 

c= |(8t,2t,6 + 5t), (A9c) 
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d = U-2 + 7t, -20 + 28t, -8 + lOt). 



(A9d) 



The maximum packing fraction is given by Eq. (11). This 
packing is vahd for 0.5777 ~ ts < t < te = |. There is no 
analytical solution for t^ because there is no analytical 
formula for (/)4,2(^)- 



5. The packing P25 



There are a total of ten intersection equations: 



E[op, oq] n (E[sr, sp] + a + b) / (/ 

£;[qo, qr] n {E[ps, pr] + b - c) / (/ 

F[o,p,q]n(d-F[o,p,q]+a)/(/ 

F[o, q, r] n (d - F[o, q, r] + b) / (/ 

F[o, r, p] n (d - F[o, r, p] + c) / (/ 

F[s, q, r] n (d - F[s, q, r] - a) / (/ 

F[s, r, p] n (d - F[s, r, p] - a - b + c) / (/ 

F[s,p,q]n(d-F[s,p,q]-b)/(/ 

F[qo, qs, qp] H (d - F[qo, qs, qp] - c]) / (/ 

F[qo, qs, qp] n (d - F[qo, qs, qp] + a + b - c) / (/ 

The lattice vectors and the offset are given by 



^(130 + t, 88 - 20t, 118 - 149t), 
^(94 - list, 136 - 92t, 106 + 37t), 



a = 
b = 

c=^(250 
d=^(-38 



179t,-44 + 10t,130 + t), 
-29t,100-94t,-62 + 65t). 



(AlOa) 
(AlOb) 
(AlOc) 
(AlOd) 
(AlOe) 
(AlOf) 
(AlOg) 
(AlOh) 
(AlOi) 
(AlOj) 



(Alia) 
(Allb) 
(Allc) 
(Alld) 



The maximum packing fraction is given by Eq. (12). This 
packing is valid for | = te < t < tr = 3 — |Vl05 ^ 
0.9506. tr is a root of the quadratic equation btj — 30^7 + 
24 = 0. 



6. The packing P26 

There are a total of nine intersection equations: 

E[po, pq] n (E[rs, rq] + a - b) / 0, (A12a) 

£;[qo, qr] H {E[ps, pr] + b - c) / 0, (A12b) 

E[rs, rp] n (£;[qo, qp] + c - a) / 0, (A12c) 

F[o, q, r] n (d - F[o, q, r] + b) / 0, (A12d) 

F[o, r, p] n (d - F[o, r, p] + c) / 0, (A12e) 

F[s, p, q] n (d - F[s, p, q] - a) / 0, (A12f) 

F[s, q, r] n (d - F[s, q, r] + b - 2a) / 0, (Al2g) 

F[op, oq, or] n (d - F[op, oq, or] + a) / 0, (A12h) 

F[sp, sq, sr] n (d - F[sp, sq, sr] + c - 2a) / 0. (A12i) 

The lattice vectors and the offset vector are given by 

a = ^(88 - 20t, 118 - 65t, 58 + 25t), (A13a) 

b = ^(-104 + loot, 46 - 5t, 106 + 25t), (A13b) 

c = ^(-56 + loot, -146 + 115t, 34 + 85t), (A13c) 

d = ^(156 - 130t, 90 - 55t, -30 + 5t). (A13d) 



The maximum packing fraction is given by Eq. (13). This 
packing is valid for 0.9506 ^ tj <t <ts = I. 

Appendix B: Packings of truncated triangular 
bipyramids with four particles per unit cell 

For the two four-particle packings P4I and P42 the 
polyhedron X as well as the lattice vectors a, b, and c are 
as defined in Appendix A. The other three particles are 
given by —X and ±MX with the transformation matrix 
M given by 



M 



1 — 2x^ —2xy —2xz 
—2yx 1 - 2y'^ —2yz 
-2zx -2zy 1 - 2z^ 



(Bl) 



Mb = b - a, 


(B2a) 


Md = f-e, 


(B2b) 


ac = 0, 


(B2c) 


af=0. 


(B2d) 



Here {x^y^z) = {ax^ay^az)/\\3i\\ is a unit vector in the 
direction of a. M describes a reflection across a plane 
that passes through the origin and is perpendicular to a. 
The offsets between X, —X, and ±MX are given by d, e, 
and f, respectively. The optimization problem is then 
solved with the intersection equations for the particular 
packing alongside the constraints 

Ma = -a, 
Mc = c, 

a • (2b - a) = 0, 
a- (e-d) =0, 

For the packing P4I there are a total of ten intersection 
equations: 

F[o,q,r]n(F[ps]+a)y^0 (B3a) 

^[oq, or] n {E[sp, sq] + b) 7^ (B3b) 

F[o,r,p]n(d-F[o,r,p])/0 (B3c) 

F[s, q, r] n (d - F[s, q, r] + a - b) ^^ (B3d) 

F[s, r, p] n (d - F[s, r, p] - b) / (B3e) 

V[qo] n (e + MF[o, q, r] + a) / (B3f) 

F[o, q, r] n (e + MV[qo] + a + c) 7^ (B3g) 

F[s, p, q] n (e + ME[to, rp]) ^ (B3h) 

£;[ro, rp] n (e + MF[s, p, q] + c) 7^ (B3i) 

^[o,q]n(f-M^[o,q])7^0 (B3j) 

The same intersection equations are satisfled in the pack- 
ing P42. Additional there are the following two intersec- 
tions: 

y [qs] n (f - MF[s, p, q] + a - b) / (B4a) 

F[s, p, q] n (f - My [qs] - b) / (B4b) 

Because of the nontrivial form of M, Eq. (7) cannot be 
solved analytically for four-particle packings. All solu- 
tions are obtained numerically. The packing P4I is valid 
for 0.5010 ^ts<t<U^ 0.5321 while P42 is valid for 
0.5321 ^t4<t<t5^ 0.5776. Note that there is no ana- 
lytical solution for ^3, t^ and t^ and there is no analytical 
formula for (/)4,i(t) and (/)4,2(^)- 
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